UNILE-CBR2 

Bounds on the Sum of a Divergent Series 



Rajesh R. ParwaniQ 



Department of Physics, University of Lecce, 
Via Arnesano, Lecce, Italy. 

September 30, 2000 
Abstract 

Given a truncated perturbation expansion of a physical quantity, 
one can, under certain circumstances, obtain lower or upper bounds 
(or both) to the sum of the full perturbation series by using the Borel 
transform and a variational conformal map. The method is illus- 
trated by applying it to various mathematical toy-models for which 
exact results are known. One of these models is used to exemplify 
how non-perturbative contributions supplement the sum of a Borel- 
nonsummable series to give the final exact and unambiguous result. 
Finally, the method is applied to some physical problems. In particu- 
lar, some speculations are made on the phase of quantum electrody- 
namics at super-high temperatures from a study of its perturbative 
free-energy density. 



1 Introduction 

After more than half a century, perturbation theory is still the best analyt- 
ical tool we have for computations in quantum field theory. Unfortunately 
in many cases the perturbation parameter, A, is not small. So the series, of 
which in practice only a few terms are known, diverges or gives a poor rep- 
resentation of the physical quantity. Consequently, several techniques have 
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been used to estimate the full sum of a series from the small number of 
given terms. In this paper, one particular method which achieves this goal 
is described in detail. The method was introduced in Ref.(l|]. 
Suppose one is given the perturbation expansion 

N 

S N (\) = J2fn^- (1) 
n=0 

Since perturbation expansions in quantum field theories have the generic 
behaviour f n ~ n! for n large [[|, it is natural to consider the Borel transform 

N f 

^)=E7' (2) 

n=0 U - 

The series (0) can then be recovered from the Borel integral, 

1 r 00 

SnW = - e^ x B N (z)dz. (3) 
A io 

If the exact function B(z) = B^z) is known, Eq.(|3]) can be taken as defining 
the exact sum, S = Sqo, of the full perturbation series if the Borel integral is 
well-defined. The poor convergence of ([!]) is then attributed to an expansion 
of B(z) in a power series as in @ and the subsequent use of that series 
beyond its radius of convergence in (|3|). 

Suppose B(z) has only one singularity in the complex z-plane (the Borel 
plane), at z = — 1/p, with p real and positive. Then the radius of convergence 
of the Borel series (0) is 1/p. Therefore in order to reconstruct an approxi- 
mation to the exact sum, one must extend the domain of convergence of the 
partial series in Eq.(0). The method of Loeffel, Le-Guillou and Zinn- Justin 
fH, f|] is to use a conformal map 

VI + zp - 1 

w{z) = — , 4 
Vl + zp + 1 

which maps the z-plane into a unit circle in the w-plane, with the singularity 
at z — —1/p mapped to w — — 1. The inverse of this map is 

z( w ) = — T 1 (5) 
p (1 — wy 

Making the change of variables (|5]) in (|3|) one obtains, 
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1 r 1 d? 

S N (X) = - dw — e- z ^' x B N (z(w)) , (6) 
A Jo aw 

where Bpf(z(w)) is understood as a power series in w obtained by an expan- 
sion of (^). In terms of the variable w, the series B^(z(w)) converges for 
\w\ < 1, so that the potential problem in (|]) is only at the upper limit of 
integration. Now, the difference between B N (z(w)) and B N+ i(z(w)) begins 
at order tu^ 1 . If one knows the coefficients f n only up to n = N, then it is 
consistent to keep only terms up to order w N in B^(z(w)). Performing this 
truncation in (|6|) and riverting back to the z variable, one obtains 



q m _ 1"^ fn (iY (2n + k-l)\ r _ Z/A {k+n) 

Sn{X) is a non-trivial resummation of the original series Sjv(A). In Ref.||, 
Eq.(|?p has been used to resum long perturbation expansions for critical 
exponents, with z — — 1/p the location of the instanton singularity of 0| field 
theory. Note that after a scaling, (0) may be written as 



n=0 



n\ \p ) r^L k\(2n — 1)! jo 



CO 



MA) = ff ^ f-V £ ( n,o ^"im 1 /" ^e-^(A,)^«). (8) 



Since u>(z) is a bounded and slowly varying function, this shows that the 
resummed expression Sjv(A) is a much slower varying function of the coupling 
than the original series Sjv(A). 

In recent years, the resummation ([?D has also been applied to quantum 
chromodynamics (QCD), with z = —1/p the location of the first ultravio- 
let renormalon pole || [7]]. Actually, QCD is Borel-nonsummable ||, which 
means that B(z) also has poles for z > 0, rendering the Borel integral ambigu- 
ous. In such a case, one can still use (|3|) to define the sum of the perturbation 
series once the integral is made definite through some prescription such as 
the principal value. 

In all the applications of (0) in |4], ^ |6|, p is a known and fixed constant 
which determines the location of the singularity of B(z) closest to the origin. 
In Ref.[|J, the expression (J7|) was used as the starting point for a new tech- 
nique which can be used even in cases when the singularity structure of B(z) 
is unknown. In the following section, I describe in detail how the technique 
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of [Tj can be used to obtain lower or upper bounds for some series, and also 
in many cases to obtain accurate estimates of the exact sum. Then in Sec. (3) 
the procedure is illustrated with several toy models. In Sec. (4) the method is 
used to study a toy model for Borel-nonsummable series, and the difference 
between the sum of the perturbation series and the exact quantity, which 
includes non-perturbative contributions, is emphasized. In Sec. (5) I discuss 
how one can simultaneously obtain both upper and lower bounds on the 
sum of certain series. Some physical applications are discussed in Sects. (6- 
8). In Sec. (6), the perturbative Euler-Heisenberg expansion for the effective 
action of quantum electrodynamics (QED) in a background magnetic field 
is resummed and the result compared with Schwinger's exact expression. In 
Sec. (7), the infrared fixed point of a <p\ theory is determined from the long se- 
ries for its beta function and the results compared with those in the literature. 
In Sec. (8) the free-energy density of super-hot quantum electrodynamics is 
studied and some speculations made. A summary and the conclusion are in 
Sect. (9). Finally, the Appendix contains an approximate but very important 
analysis of the relevant equations which explains the empirically observed 
trends. 

2 The Variational Conformal Map 

The conventional transformation of a series 

N 



S N (\) = J2fn^- (9) 
n=0 

by the Borel-conformal map method results in the reorganised expression 

(10) 

with p a fixed constant. However, notice that p does not figure in (|[) but 
enters (ITU]) only through the conformal map (HI). Thus instead of treating p 



as some fixed constant as in Ref. ||, [|, ||, [7|, one may consider p > a free 
parameter which defines the conformal map (H) and Eq.(^DJ) then represents 
a continuous family of resummations, one for each value of p. Thus, from 
now on, p will not refer to the location of some singularity in B(z). In fact, 
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absolutely no knowledge about the singularity structure of B(z) is required 
for the method elaborated below. 

Of course having liberated ourselves from the usual interpretation of p 
in (|H]), we also lose defmiteness in our resummation. Therefore some new 
condition must be imposed to fix the value of p and hence of the expression 
(|T0|). For a start, for each N, choose p > to be the location of an extremum 
of Sn(\,p) flTJ. Since (|T0|) depends on A, the value of p in principle will also 
depend on A. However the procedure then becomes too unwieldy, and to 
simplify it p is determined at some reference value A = Ao, say at the mid- 
point of the range of interest: 

dS N (X ,p) _ Q 
dp 

As mentioned after (H), S(X,p) is a relatively slowly varying function of A, 
hence determining p at some fixed A = Ao, and then using the same p in (|lOl ) 
for various A is sufficient for practical purposes. Indeed, in many applications, 
one requires the sum of the series at one particular value of the coupling, or 
in a very narrow range, so the simplification made in ([IT]), is both useful and 
sufficient. 

In some cases, (^JJ) will not have any solutions. For example, if all the /„ 
are of the same sign, Sn(X,p) will be a monotonic function of p. Thus for 
the procedure to work, at least some of the f n must be of a different sign. 
This will be assumed to be the case from now on. Suppose next that f\ ^ 0. 
Then for p large, and since the p dependence of w(z) is mild, 

S N (p)~fo + j- (12) 

Therefore for /i < 0, Sn will first decrease as 1/p increases and then increase 
when the next term f m /p m > dominates the sum. If the last non-zero 
term Jn/p N is positive then one expects Sn(p) to have a global minimun at 
some p > 0. However if fN/p N < then Sn(p) can become very small as 
p — > + and so the minimum would likely be only a local extremum. From 
this heuristic argument one concludes that if the first non-zero coefficient 
f n , (n > 0) is negative, Eq.(|TT|) will have global minima solutions for some 
N. 

For conciseness, unless otherwise stated, I will discuss from now on only 
the case when a global minimum exists, as the arguments for the global 
maximum case are then obvious. For a series Sn(\,p), define p(N) to be 
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the position of the global minima for the values of N when they exist, and for 
other values of N let p(N) denote the location of the local minima. Also, let 

S N = S N (\,p(N)). (13) 

Obviously the definition of p(N) has been chosen because it is useful. If 
p(N) is a location of a global minimum, then Sn certainly is a lower bound on 
Sn{^,p)- However this does not imply that Sn is a lower bound on the exact 
sum S of the pertutbation series because for finite N, S might lie outside of 
the space of resummations labelled by p. For each N, let p*(N) be the value 
of p that is optimal, that is, it is the value which when used in (|10D gives the 
best estimate of S. Define, S N = Sn{^,P*{N)). Then for a global minima 
one has 

S N < S* N (14) 

Presumably S N converges to S as N — > oo. Then for those N when global 
minima exist, 

S N ^oo<S. (15) 

(This implicitly assumes that the sub-sequence of global minima is infinite: 
That is, given any positive integer N , there is some n > N for which f n is 
positive.) 

Though the inclusion of global minima in the definition of p(N) is quite 
clear, the inclusion of local minima requires some explanation. As mentioned 
earlier, if Sn(^-p) has a global minimum, then Sn+i(\p) will probably not 
have a global minimum if /jv+i is negative because Sn+i(\p) might become 
very small as p — > + . However Sn+i{^,p) will still have a local minimum 
(and hence also a local maximum). The local minimum will occur for mod- 
erate values of p close to p(N), the global minimum position of Sjv(A,p). 
Therefore one might expect that the local minimum for Sn+i{^,p) still gives 
a bound on S N+1 . Indeed, suppose that the inequality 

Sn < S N +i (16) 
holds for all N up to N = oo, then clearly 

S N < S (17) 

and one concludes that Sn(\,p(N)) is a lower bound on the sum of the 
full perturbation series. Of course proving ([16]) is equivalent to knowing f n 
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explicitly for all n, which is not the situation in reality. As a practical matter, 
it is sufficient to draw a conclusion after observing the trend (|T6| ) (or lack 
of) for the available terms of the series. Arguments given in the Appendix 
indicate that the trend, once started, will continue. 

An interesting question is whether the inequality (|15D is saturated. In 
most of the examples studied this has been found to be the case. In fact 
the convergence is so rapid that one can conclude with a high degree of 
confidence, not only that Sn is a lower bound, but that it is close to the 
exact value. However a toy model in Sec. (3.6) and a physical example in 
Sec. (7) show that the series Sn, though forming a bound on the actual value 
S, might not converge to S. In both these examples the exact result S(X) 
had a first derivative dS(X)/d\ that varied rapidly with A. An explanation of 
why in those cases the bounds Sn do not converge to the exact value is given 
in the Appendix. However, as described below, such situations can also be 



taken care of by the resummation procedure (|T0| - |TT 



Recall that when a local minima occurs for some N, one expects also a 
local maximum. Define Pn to be the position of those local maxima and Sn 
the corresponding value of the resummed series. It turns out that Sn also 
satisfies an inequality like However since for the sequence Sn of only 

local extrema one does not seem to have a statement like ( |I1| - |15| ), it is not 
a priori obvious that they form bounds to the exact result. Of course if one 
believes that the p*(N) as defined above always take moderate values, then a 
statement like ( |HD might also be made for the sequence Sn- Furthermore, in 
Sec. (5) I describe how an auxiliary series S N (X,p) defined from Sn(X,p) can 
be used to obtain upper bounds to S, when Sn gives lower bounds. Using 
the sequences Sn, Sn, and S N , one can obtain constraints on S and, with 
some physical input, an estimate of S itself. In the Appendix it is explained 
why the alternate bounds formed from Sn and the auxiliary series give better 
approximations to the exact value S itself when dS/dX varies rapidly with 
A. 

Although the main discussion in this paper will be for the p(N) (or p(N) 
andp'(N)) as defined above, in some cases one finds (by inspection) solutions 
Po(N) to ( |ITD which have the property that as iV — > oo, Po(N) — > p , a 
constant. In such a case one is tempted to speculate that the fixed point po 
indicates the existence of a singularity at z = — 1/p in B(z). This is indeed 
found to be the case in the examples studied, though apparently the converse 
is not necessarily true: the singularities of B(z) need not show up as solutions 
of (fnp . Furthermore, the convergence of the series S^n is not expected to 
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be monotonic since in general the po(N) refer to both maxima and minima. 

Though not manifest at first sight, the analysis in the Appendix suggests 
that even the p(N) defined above Eq.fllBD actually will converge to a fixed 
value as N —>■ oo. This trend can be observed in the examples studied. 



3 Mathematical Models 

In this section a number of mathematical models are studied to illustrate 



the resummation technique of Eqs. (|10| - |iTD . The models all define Borel- 
summable series, that is, B(z) has no singularities at real, positive z. A Borel- 
nonsummable example will be considered in Sec. (4). All of the examples 
here have global minima as solutions to the extremum condition fllTD for 
some N. (From any such series Sn one can trivially construct C — Sjy, with 
C a constant, which then gives an example of a series with global maxima 
solutions to (]TlD ). The reader is reminded that the existence of global minima 
does not by itself imply that one has obtained a lower bound on the sum of 
the series. The additional condition ([TJ]) must be satisfied for a lower bound. 
Example (5) below shows how one can end up with an upper bound from 
global minima ! (See also the Appendix). 

It is also re-emphasized that although in these toy models the exact sin- 
gularity structure of B(z) is known, that information will not be used in the 
resummation. The resummation of the partial series 

N 
n=0 

will proceed using Eqs.(|TU]-|n|). The exact B(z) will only be used to compare 
the resummed results with the exact sum of the full series given by the Borel 
integral 

1 r°° 

S(\)= e - z/x B(z)dz. (19) 



A Jo 



3.1 B(z) = ^ 

The first model is defined by the Borel function, 



B{z) = T L-. (20) 
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Expanding B(z) to iV-th order in z and using the result in (Q) gives the 
truncated series 



N 

^ = E(- A ) n ^ ! - (21) 

n=0 

The divergent nature of this series is displayed in Fig. (la). Using f n = 
(— l) n n\ in ( |TUD and solving Eq.([TTp at the reference value Ao = 1 gives 
the following solutions (minima): p{2) = 2.65, p(3) = 5.1 and p(4) = 8.4. 
As expected from the arguments of Sec. (2), there is no solution for N = 
1, the solutions for N = 2 and N = 4 are global minima while that for 
N = 3 is a local minimum. The resummed series is shown in Fig. (lb). The 
convergence of the Sn is monotonic and satisfies the condition flip]). Hence 
the approximants Sn can be argued to form lower bounds to the exact result. 
That this is indeed the case can be seen by inspection of the exact result in 



Fig. (lb) as obtained from ( 20|) and (|T9|) . Furthermore, it is clear that the 
lower bounds converge to the exact value, and so may be used to estimate 
it. At A = 0.5, the exact value is 0.722657, while the approximants are 
S 2 = 0.704, S 3 = 0.709, S 4 = 0.711. 

For this model there are other solutions, for each N, to the extremum 
condition (|ll]) in addition to the minima. By inspection one picks out the 
sequence of values, Po(3) = 1.6 (local maximum), po(4) = 1-3 (local mini- 
mum) and Po(5) = 1.15 (local maximum) as plausibly approaching a fixed 
point. Indeed we already know that in this model the exact singularity of 
B(z) is a pole at z = —1/p— —1, so the fixed point probably refers to the 
location of this singularity. It should come as no surprise that if the approx- 
imants ([LC]) are evaluated at the values p (N), the convergence to the exact 
value will be much faster, and this is indicated in Fig.(lc). At A = 0.5 the 
values of the approximants are S^ 3 = 0.726, S , ( )4 = 0.7219, = 0.7228. 
Notice however that in this case the convergence is not monotonic (and was 
not expected to be). 

The alternate sequence Sn for this model will be discussed in Sect. (5). 

3.2 B(z) = ^ 

The second model is defined by the Borel function, 



«w = rh (22) 
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which has poles only on the imaginary z-axis. The truncated perturbation 
series corresponding to (|22|) is 

N 

&r=£(-l) n X " n (2n)!. (23) 

n=0 

The divergent nature of this series is displayed in Fig. (2a). Using f n = 
(— l) n (2n)! in (|T^) and solving Eq.(|Il]) at the reference value Ao = 1, gives 
the following solutions (minima): p(4) = 1.45, p(5) = 2.9, p(6) = 4.5, p(7) = 
6.3, p(8) = 8.4. Notice that although /Wh-i = 0, the approximants S 2 n+i do 
exist and are different from SW' This is because of the way the resummation 
is performed in Eq . (fl0|) . As expected on general grounds, the minima for 
N = 4, 5, 8 are global minima while those for N = 6, 7 are local minima. 

The resummed series is shown in Fig. (2b) together with the exact result. 
The convergence is again monotonic, satisfies the condition (|T6"D, and the 
approximants Sn, are indeed lower bounds to the exact result. Furthermore, 
compared to the divergent series (P3"D, the resummed series is close to the 
exact result for moderate values (~ 0.3) of the coupling. At A = 0.3, the 
exact result is 0.89, while the approximants give, S 4 = 0.86, 5*5 = 0.868, S 6 = 
0.871, S 7 = 0.873, S 8 = 0.874. 

3.3 B(z) = eT z 

The third model is defined by the Borel function, 

B(z) = e~ z . (24) 

This function is regular everywhere in the Borel plane and the exact sum of 
the perturbation series has the closed form S(\) = 1/(1 + A). Thus in this 
case a power expansion of S(\) is actually convergent for |A| < 1. However 
close to A = 1, the convergence is very slow and one requires a large number of 
terms of the series to obtain accurate results. The utility of the resummation 
procedure even in this case is now demonstrated. 

The truncated perturbation series corresponding to ( p4|) is 

N 

SW = £(-l) n A™. (25) 
o 

Notice that there is no factorial growth of the coefficients in (|25|). The slow 
convergence of this series for A close to 1 is displayed in Fig. (3a). Using 
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f n = (— l) n in (|TUD and solving Eq.(|TTD at the reference value A = 1, gives 
the following solutions (minima): p(2) = 1.2, p(3) = 2.9, p(4) = 5. The 
resummed series is shown in Fig. (3b) together with the exact result. The 
convergence is rapid, monotonic, satisfies the condition ([16]), and the approx- 
imants Sn form lower bounds to the exact result. With only four terms, the 
resummed series already shows a dramatically improved convergence com- 
pared to the original series (|25|) , even for couplings as large as A = 0.8: 
The exact sum at that coupling is 0.556, while the resummed values are 
S 2 = 0.492, S 3 = 0.507, S 4 = 0.512. 

This example illustrates that the resummation procedure ( |TQ] - prT| ) is useful 
even for a convergent series, especially when one is close to the radius of 
convergence and when only a few terms of the series are available (which is 
often the situation in practice). Furthermore this example, and the last one 
in Sec. (3.2), emphasize that the parameter p has no obvious relation to the 
singularities of B(z). 

3.4 B(z) = -^ + -h 

The fourth model is defined by the Borel function, 

B ( z ) = TT ^ ■ ( 26 ) 

l + z 1 + 



The truncated perturbation series corresponding to (|26|) is 



N 

Sn = T, f n \ n , (27) 







with 



f 2n+1 = -(2n+l)!, (28) 
f 2n = (l + (-l) n ) (2n)L (29) 

This is is an example of a series which is not strictly alternating. The diver- 
gent series is displayed in Fig. (4a). Using the values of /„ given above and 
solving Eq. (|TTD at the reference value Ao = 1, gives the following solutions 
(minima): p(4) = 1.5, p(5) = 2.8, p(Q) = 4.5, p(7) = 6.4. As expected from 
the signs of the f n , there are no solutions for N < 3, the solution for N = 4 
is a global minimum, while those for iV = 5,6,7 are local minima. Again 
note that exists and does not equal S 5 eventhough f 6 = 0. 
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The resummed series is shown in Fig. (4b) together with the exact result. 
The convergence is monotonic, satisfies the condition (|16D, and the approxi- 
mants Sn form lower bounds to the exact result as expected. The resummed 
series shows a much improved convergence compared to the original series 
(|27|), even for couplings as large as A = 0.5: The exact value is 1.521, while 
the resummed approximants give, S4 = 1.397, S 5 = 1.426, 5*6 = 1.437, S 7 = 
1.443. 

In this example, no sign of the sequence Po{N) was found even though 
there is a pole of the Borel function B(z) at p = —1. This is taken as support 
of the statement at the end of Sec. (2). Actually, a comparison of the p(N) 
values of this model with those of the model in Sec. (3.2) shows that they are 
almost identical ! This means that the 1/(1 + z 2 ) part of the Borel function 



in (|26|) is dominating the behaviour of p(N) and thus masking the singularity 
at z — — 1. 

This example also displays the following curiosity: Although the series 



(27) is divergent, if one keeps only the first two terms /o an d fi, then the 



unresummed expression S'i(A) agrees quite well with the exact result S(X) 
(See Fig. (4a)) for a large range of couplings! A similar 'accident' occurs for 
the free energy density of hot SU(3) gauge theory, where it is found that 
the usual second order perturbative contribution already agrees with the full 
lattice result (see H and references therein). 

The fifth model is defined by the Borel function, 

B{z) = , (30) 

which has a branch cut beginning at z = — 1. The truncated perturbation 
series corresponding to (BO) is 




<SW = £ hr 2n(2n-3)!! (31) 



with (—1)!! = 1. The divergent series is displayed in Fig. (5a). Solving Eq.flTTD 
at the reference value Ao = 1, gives (global minima): p(2) = 1.3, p(A) = 
1.1, p(6) = 1. No solutions were found for iV = 3 or 5. This example 
illustrates that although the heuristic argument of the Sec. (2) suggests a 
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solution to flTT|) in general, such a solution might fail to exist for some values 
of N. 

This example has a number of other peculiarities. First notice that the 
global minima sequence p{2) = 1.3, p(4) = 1.1, p(6) = 1 is identical to 
the sequence Po(N) (see the end of Sec. (2)) which appears to converge to 
the fixed point po = 1, and which happens to locate the exact singularity of 
B(z). 

The resummed series is shown in Fig. (5b) together with the exact result. 
The convergence is rapid and monotonic but the approximants Sn do not 
form lower bounds to the exact result as might have been naively expected 
because S 2 > 5* 4 > S e , which is opposite to what is required. Indeed, as the 
figure shows, the approximants approach the exact result from above and 
hence appear to form upper bounds eventhough the p(N) are positions of 
global minima! An explanation of this oddity is given in the Appendix. See 
also the discussion prior to Eq. (|14|) . 

For A = 2, the exact result is —1, while the succesive approximants give, 
S 2 = —0.907, 5*4 = —0.986, S 6 = —0.996. Thus the approximants converge 
to the exact value. 

3.6 S{\) = -sin(TrA) 

The sixth model is defined by the exact expression 

S{\) = -sin(vrA) . (32) 
The partial series corresponding to (|32|) is 

This series is convergent but it has been chosen to test the resummation pro- 
cedure for cases when the exact expression S(X) is not a monotonic function 
of the coupling. The partial series is displayed in Fig. (6a). 

Solving Eq.([0]) at the reference value Ao = 1, gives (minima): p(3) = 
1.5, p(4) = 3.1, p(5) = 5, p(6) = 7.3, p(7) = 10. Notice that solutions exist 
for N = 4 and N = 6 (and are different from those for N = 3, 5) eventhough 
f^ — fe — 0. (See the discussion for the example in Sec. (3. 2).) 

The resummed series is shown in Fig. (6b). Though the convergence of the 
S*7v is rapid and monotonic, and although the approximants Sn form lower 
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bounds to the exact result, they do not give a good estimate of the exact 
result when A > 0.2. This seems to be because the resummed expressions 
Sn are monotonic functions of A while the exact expression is not (see the 
Appendix). 

Now, recalling the general discussion of Sect. (2), one expects there to be 
local maxima solutions to Eq.([0]) when N = 5 and N = Q. Indeed the 
solutions are p(5) = 0.15 and p(6) = 0.275. The curves for the approximants 

and Sq are shown in Fig. (6c). These are seen to form upper bounds to 
the exact result and also give very good estimates ! Thus it seems that 
one should use all the solutions to (|i~T|), that is the p(N) and p(N), to get 
all possible bounds on the sum of the series. Then using some additional 
physical or theoretical information or prejudice, one can decide near which 
of the bounds (upper or lower), the exact result lies. A physical example is 
given in Sec. (7). 

In this example, eventhough we used a series up to = 7, only for two 
values could we form the upper bounds. In Sec. (5) I describe how one can 
get additional upper bounds from the series by first constructing an auxiliary 
series. 



4 A Borel-Nonsummable Model 

For most of the physical quantities calculable from the Standard Model of 
particle physics, the perturbation expansion is not expected to be Borel 
summable. In simple terms, this means that the function B(z) has poles 
on the positive semi-axis of the Borel plane thus rendering the Borel integral 
(P~9| ) ambiguous. If one choses an ie presciption then the resulting ambigu- 
ity is in the imaginary part. Sometimes these imaginary parts are of direct 
physical relevance || 1Q| . More generally they indicate that the perturbation 



series does not give the full answer, but must be supplemented with some 
non-perturbative contributions || ||. Since the imaginary parts will be of 
the form e _1//A , the additional real non-perturbative terms are expected to 
take the same form. These expectations have been confirmed in some lower 
dimensional field theories (see || and references therein). 

One can also construct mathematical models to illustrate the arguments 
of the last paragraph. Suppose, for simplicity, that the only singularity of 
B(z) for z > is a single pole at z — q. Then the sum of the perturbation 
series can be defined by the principal value prescription, 
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S pert = X V dz e- z/A B(z) . (34) 
A jo 

With this definition, one focuses only on the real part of the physical quantity. 
Suppose furthermore, again for simplicity, that B(z) is integrable at infinity. 
Then the Eq.(|34T) may be rewritten as 



Spert = -r / dz {e- z ' x - e-«/ A ) B{z) + —-V \ dz B(z) . (35) 
A Jo A Jo 

The first term on the right-hand-side is finite and unambiguous. Call it S exact , 
and denote the second term on the right-hand-side as S np . Thus we have 

Sexact(A) = S per t(\) — S np (\) . (36) 

In this way, the exact result S exact , has been broken into two components, 
Sp er t which is purely perturbative, and S np which is purely non-perturbative. 
However while S exact is well-defined, both S pert and S np are only defined 
through the principal value prescription. Though highly simplified, this 
model plausibly represents the situation, for example, in Quantum Chro- 
modynamics (QCD). 

In order to illustrate explicitly the resummation technique ( |T0| - |Tl~| ) in the 
Borel-nonsummable case, set 

B & = n ■ L t • ( 37 ) 

(1 + z)(5 - z) 

Then, 

SpertW = -V / ~ ; r~ r . (38) 



A Jo (l + z)(B-z 
The truncated perturbation series corresponding to (|38D is 

N 

S N = Y,fnX n , (39) 



with 

u = ((-ir+ ^ n+X) ) ^. (40) 

o 

The divergent series is displayed in Fig. (7a). The solution of Eq.([TT|) at the 
reference value Aq = 1, gives (minima): p(2) = 2.8, p(3) = 5.4, p(4) = 
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9, p(5) = 13.5. The resummed partial series is shown in Fig. (7b) together 
with the exact sum of the full perturbation series given by fl3~8|) . The conver- 
gence is rapid, monotonic, satisfies the condition ([16]), and the approximants 
Sn do indeed form lower bounds. Also, the convergence is so fast that, except 
for a small interval at intermediate coupling, S§ is not only a lower bound, 
but also a very good approximation to S pert . 

Now, corresponding to the Borel function ([37D , one has from the defini- 
tions before Eq. (|36|) , 

S np W = ^ e- 5 ^ (41) 

and 

A jo [l + z){b — z) 

In Fig. (7c), the curves for S exact and S 5 are plotted for a bigger range of 
A. This shows that except for a small range of couplings, the resummed 
perturbative approximant S§ lies above and deviates significantly from the 
exact result although, as seen above, S5 does form a converging lower bound 
to the perturbative component S per t of the exact result. The difference is 
of course the non-perturbative piece S np . In the same figure there is a plot 
of (S 5 — S np ) which, according to the definition (|36"D should approximate 
S exact . Indeed the agreement is very good and, amusingly, it improves at 
large coupling: At A = 5, S exact = 0.0533 and S 5 — S np = 0.0457, while at 
A = 10, S exact = 0.0219 and S 5 - S np = 0.0193. 

This toy model discussion illustrates concretely the arguments given in 
Ref. for the free-energy density of thermal 577(3) gauge theory. There it 
was found that the exact result, given by lattice data, differed significantly 
from the resummed perturbative result, and it was argued that the differ- 
ence was caused by the Borel-nonsummability of the theory. Assuming a 
non-perturbative component of the form y e~ Q ^ x , the constants A and q 
were determined from the difference between S exact and S 5 . A more detailed 
discussion of the results in U and their extension to full QCD and other 
non-Abelian gauge theories will be presented in pT|| . 

This simple toy model exhibits the following curiosity. In Fig.(7d) the 
second order resummed perturbative result 5*2 is plotted together with the 
exact expression S exact . The curves agree very well over a wide range of 
couplings! This shows that sometimes a resummation of purely perturbative 
results can accidently give agreement with the exact value which contains 
both perturbative and non-perturbative components. 
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5 Lower And Upper Bounds 



Most of the discussion so far has been for the approximants Sn formed from 
the p(iV)'s which are positions of global or local minima. However as men- 
tioned near the end of Sect. (2), when local minima exist, so generically must 
local maxima. Those local maxima, located at p(N), form the series S(N) 
which provides additional information. In the S(X) = — sin(A7r) example 
studied earlier, it turned out that though the Sn and Sn approximants sep- 
arately formed converging lower and upper bounds, the exact result was 
closer to the upper bounds. In the next subsection, another example, that 
of Sect. (3.4), is re-investigated to obtain upper bounds from the positions of 
its local maxima. 

Then in the second subsection, another technique, that of using an aux- 
iliary series, is introduced so that alternative bounds to a series may be ob- 
tained also from positions of global maxima. The reason why this alternative 
method is important is that, as explained in Sec. (2), a priori the sequence 
Sn does not seem to obey a condition like (|I^ , |T5[) . Therefore in practical 
problems where the exact result is unknown, any conclusions drawn on the 
basis of Sn should be confirmed by other means, such as the auxiliary series 
method. 

As explained in the Appendix, the bounds formed from the Sn are gener- 
ically monotonic functions of A and therefore do not approximate the exact 
result 5(A) very well if the latter is not monotonic. In that case, the bounds 
formed from the local extrema, Sn, and from the auxiliary series, are ex- 
pected to give better approximations to S (see Appendix). 

5.1 Bounds from Local Maxima 

Generally, for a short series which mainly has global and local minima so- 
lutions to (|ll"D , the number of values for which local maxima solutions exist 
will be even less. However for a series which is not strictly alternating, there 
will be more local extrema solutions (see the general discussion in Sec. (2)), 
and hence an opportunity to form a longer series Sn- An example of this is 
provided by the toy-model of Sect. (3.4), with exact Borel transform 




(43) 



The perturbation series corresponding to ( f43"D is 
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Sn = J2 , (44) 



with 

f 2n+ i = -(2n+l)! (45) 
/an = (l + (-l) n ) (2n)!. (46) 

From the arguments of Sect. (2), one expects local maxima to exist for 
N = 5,6, 7. An explicit analysis confirms this and gives their location. Using 
the reference point A = 1, the solutions to fllTD are: p(5) = 0.3, p(6) = 
0.7, p(7) = 1.5. The corresponding curves for Sn are shown in Fig. (8) and 
clearly form upper bounds to the exact result. At A = 0.5 the values are 
Ss = 1.97, Sq = 1.80, Sj = 1.68, while the exact value is 1.521. Note the 
expected large curvature of the bounds. 

5.2 Bounds from an Auxiliary Series 

Given a series which has global minima solutions to (|i~l"D, then by 'removing' 
the first nontrivial coefficient of that series one forms an auxiliary series which 
has global maxima solutions. In this way one can form alternative bounds 
to the sum of the original series complementing those obtained from Sn- 

To illustrate this concretely, re-consider the toy example Sec (3.1) with 
B(z) = 1/(1 + 2). The original truncated series which leads to global minima 
and lower bounds is 

N 

&r = E(-A)"n!. (47) 

n=0 

Define the auxiliary series 

S' N = {Sn ~ 1} = £ (-1)" +1 {n + 1)! A n . (48) 

A n=0 

Since the auxiliary series S' N begins with a positive non-trivial coefficient, 
it will give maxima solutions to flTl"|). The solutions at the reference point 
A = 1 are, p'(2) = 3.9, p'(3) = 8, p'(A) = 13.25, j/(5) = 20, p'(6) = 28. 
The curves for 1 + XS N (X) are plotted in Fig. (9). Indeed they are seen to 
form upper bounds to the exact result. However the convergence to the exact 
result is slower than that of the lower bounds Sn considered earlier. 
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For a second example, re-consider the example of Sec. (3.6) with S(X) = 
— sin(A7r). The partial series is given by 

ff "=£fe ( - 1, " B - (49) 

Define the auxiliary series 

A ^ A (2n + l)! 

The auxiliary series begins with a positive non-trivial coefficient and so one 
will obtain maxima as solutions to the extremum condition (|iT|). At the 
reference value A = 1, the solutions are: p'(4) = 0.3, p'(5) = 0.575, p'(6) = 
0.92. The corresponding curves for XS N are shown in Fig. (10). Not only do 
they form converging upper bounds, but in this case they approximate the 
exact function very well compared to the approximants Sn used earlier in 
Sect. (3.6). 

A question now arises for practical problems where exact results are not 
known. Is it possible to determine whether the unknown exact result lies 
closer to the upper or lower bound ? Empirically, it appears that the series Sn 
converges to the exact result S(X) if d 2 S(X)/dX 2 is very small in magnitude 
in the entire range of interest, < A < Ao- Otherwise the exact result seems 
to be better approximated by the series Sn or by the auxiliary series method 
discussed above. An explanation of why the bounds Sn have slowly varying 
first derivatives, and why the bounds Sn and those from the auxiliary series 
have faster varying first derivatives is given in the Appendix. 



6 The Effective Action of QED in a Magnetic 
field 



The first physics example is Schwinger's [12) effective Lagrangian (density) 
for QED in a uniform magnetic field B, 



elB 2 



L(X) 



■sir* Jo 



00 ds / _ . . , 1 
Coth(s) - - 



(51) 



where A = e 2 .B 2 /m i , e r is the renormalized electron coupling, and m the 
electron mass. 
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Define the dimensionless quantity 



5(A) = 100 r % footh(a) - - - ~) e- st ^ (52) 
Jo s \ s 6J 



which is related in an obvious way to (pi). An expansion of ( p2|) is given by 



S{\) = 1600 J2 /nA n (53) 

?1=1 

with 

/n ~ 2n (2n + 1) (2n + 2) ' 1 j 

and where i^2 n are the Bernoulli numbers which alternate in sign, and diverge 



factorially for large n. Equation is the Euler-Heisenberg |13| expansion 
which is equivalent to a sum of an infinite number of Feynman diagrams 
consisting of one closed fermion loop with an even number of external pho- 
tons. The Euler-Heisenberg series diverges for large values of A as shown in 
Fig. (11a). Consider now a resummation of the divergent series using ( |I0~| - 
ITJ). At the reference value A = 10, the solutions to (|ll|) are (global minima) 
p(2) = 0.77, p(3) = 1.4, p(6) = 2.25. No solutions were found for N = 3,5, 7. 

The exact result of Schwinger, given by (|52"D , is plotted in Fig. (lib), 
together with the approximants S 2 , S'4, S 6 . As in the toy model of Sec. (3.3), 
the approximants form upper bounds although the p(N) are positions of 
global minima. See again the caution prior to Eq.([rj) and the Appendix. 
The convergence of the approximants Sn to the exact expression is manifest. 
For example, at A = 10, where the Euler-Heisenberg series is badly divergent, 
the Schwinger's exact value is 5'(10) = —8.056 while the estimates are, S 2 = 
-5.9, S 4 = -6.9, S 6 = -7.3. 

The Euler-Heisenberg series ( |53"D is Borel summable [0, and therefore 
provides an example of a Borel-summable series in a theory (QED) which is 
generally considered Borel- nonsummable. However it should be noted that 
the Euler-Heisenberg series (^) represents a one-fermion-loop result whereas 
the "renormalon" singularities || which signal Borel- nonsummability are ex- 
pected to arise when higher-loop diagrams are considered. That is, Borel- 
nonsummability might manifest itself when multi-loop corrections to (0) are 
taken into account. 
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7 The 0| Field Theory 



The three-dimensional 0(J\f) symmetric <p A field theory can be used as an 
effective theory (see j|, ||| and references therein) to describe the critical 
behaviour of many physical systems near a second-order phase transition. 
For this purpose, the renormalization group functions of this theory have 
been calculated to very high order. A detailed list of references can be found 
in H For example, the beta function for the polymer case, M = 0, is 
given by 

/5(A) = -A+A 2 -0.439815A 3 +0.389923A 4 -0.447316A 5 +0.633855A 6 -1.03493A 7 . 

(55) 

where A is a dimensionless coupling (usually denoted as g or g in the literature 
||). The objective is to find the non-trivial infrared fixed point, A*, of the 
theory, which is given by the zero of the beta function, 



/3(A*) = 0, (56) 
dA 1 



A* 



= uo > . (57) 



The expression ( |55| ) is divergent for large A ~ 1 where a nontrivial zero 
is expected. Using the resummation (|l(J) with Sn denoting approximants 
to the beta function, and choosing the reference point Ao = 1, one gets as 
solutions to Eq . (|TT|) (minima), p{2) = 1.3, p(3) = 3.2, p(4) = 5.6, p(5) = 
8.6, p(6) = 12.25, p(7) = 16.5. The curves are shown in Fig. (12a). They 
form rapidly converging lower bounds but intercept the A-axis only at the 
origin, thus giving only a trivial zero to the beta function. The situation is 
similar to the sin(7rA) toy model studied earlier. 

Now, since for N = 3, 5, 7 the minima are only local, one expects lo- 
cal maxima to exist also. Indeed the local maxima are located at, p(3) = 
0.18, p(5) = 0.21, p{7) = 0.19. The curves for Sn are shown in Fig(12b). 
Though these are local maxima, they obey the inequality Sn+i > Sn and 
so appear to form lower bounds ! Since the lower bounds due to Sn are 
higher than the lower bounds due to Sn (Fig. (12a)), one may argue that 
those due to Sn provide more accurate information. The curves in Fig. (12b) 
do in fact have a non-trivial zero. For N = 7, the zero is near A = 1.425. 
In this respect, the physical example here is different from the toy- model of 
Sec. (3.6): here both the approximants Sn and S(N) provide lower bounds, 
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so it is simply a matter of choosing the highest lower bound to estimate the 
sum of the series. 

Given the physical importance of this example, it is useful to perform 
further checks. So let us re-analyse the problem using the auxiliary series 
method of Sect. (5). That is, let S N (X,p) = Sn(X,p)/X. The auxiliary series 
will then obviously have maxima as solutions to ([□]). The solutions at the 
reference point A = 1 are p'{2) = 0.6, p'(3) = 0.875, p'(A) = 0.33, j/(5) = 
0.4, p'(6) = 0.2525. (Actually, p'(5) is a point of inflexion). The curves 
for XS N are shown in Fig. (12c). Again, though determined by positions of 
maxima, the curves appear to form rapidly converging lower bounds which 
have a non-trivial zero. 

In order to compare with the results in the literature, more precise num- 
bers are now quoted. Conjecturing that the curves in Fig. (12c) are indeed 
lower bounds to the exact result, the N = Q curve, shown magnified near 
its nontrivial zero in Fig.(12d), gives an upper bound of 1.4193 to the non- 
trivial zero of the beta function. Re-optimizing the N = 6 equation ( [11]) at 
A = 1.419 does not change the curve for XS N significantly to modify that 
bound at the level of accuracy quoted. The slope of the beta function at the 
non-trivial zero can also be determined from the N = 6 curve in Fig.(12d). 
It is uj = 0.7955. Since the curves appear to get steeper as N increases, this 
value of uj is a lower bound. These values can now be compared with those 
of Ref.0: there it was found that A* = 1.413 ±0.006 and uj = 0.812 ±0.016. 
The agreement with the bounds found here is excellent. (Comparison of 
results obtained by other methods and other authors may be found in (0). 

It is important to note the significant conceptual difference between the 
methodology used in this paper and that employed in Ref.p. In [|] the 
authors used a Borel-Leroy transformation with a variable parameter b but 
they used the conventional conformal map M with the value of p fixed at the 
precise location of the instanton singularity, p* = 0.166246 (the value for M = 
0). The parameter b, together with some other parameters introduced in || 
were used to check the convergence of the series and to estimate their errors. 
Here instead the usual Borel tranform (Q) is used but the conformal map 
has a single variational parameter p determined according to the condition 
(Trip. From the numbers quoted above, it is clear that the values p(N) used 
here are not the same as the value p*. Furthermore, in the approach of this 
paper, no assumption about the analyticity structure of the Borel transform is 
made except that in order for the resummed perturbation series to faithfully 
represent the physical quantity, the series should be Borel summable (which 
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I simply take to mean that there are no poles on the positive Borel axis, and 
the Borel integral converges at its upper limit). Borel summability of the <j)\ 



theory has been established in [15 



From this example it is clear that the novel resummation presented here, 
with the parameter p determined from ([LTD , can be used to complement the 
analysis done in pf. 



8 Hot Quantum Electrodynamics 

The fine structure constant a of QED is so small (~ 1/137) that in prac- 
ticethe perturbation series converges even without any resummation. How- 
ever renormalization group arguments indicate that a will increase at high 
energies. Since the growth is only logarithmic, the energy scale must be 
exponentially high before the perturbation series fails to converge. 

The free-energy density of QED at very high temperature (T) has been 
computed up to order a 5 / 2 |T^, [T7|]. The temperature was assumed to be 
high enough so that the electron mass could be neglected. Let A = (a/ix) 1 / 2 . 
Then the normalised free-energy density of QED at the MS renormalisation 
scale p, = 2tiT, is given by |L6|, [L7| 

F/F = 1 - 1.13636A 2 + 2.09946A 3 + 0.488875A 4 - 6.34112A 5 . (58) 

where F = ll7r 2 T 4 /180 is the free-energy density of a non- interacting plasma. 
Figure (13a) shows the plot of (|58"D at different orders. At large coupling 
(super-high temperatures) the series diverges, exhibiting a behaviour simi- 
lar to that of Yang-Mills theory at low-temperatures. The convergence at 
large coupling can be improved by using the resummation technique ([T0|-pT 



Using the coefficients from (0), the solutions of (|Tl"| ) at the reference value 
Ao = 0.5 are (minima): p(3) = 0.7, p(4) = 1.75, p(5) = 3. 

The resummed approximants are shown in Fig. (8b). The convergence is 
clearly much better than in Fig. (8a), suggesting that the curves not only form 
lower bounds but also good estimates to the perturbative free-energy density. 
The emphasis on 'perturbative' is because QED, like QCD, is probably Borel 
non-summable, and so the large coupling perturbative results, even when re- 
summed, might differ significantly from the exact result by non-perturbative 
contributions of the form e~ q l x j\. 

If one assumes that the potential non-perturbative contributions lower 
the perturbative result, as happens in QCD (see [l[] and references therein) 
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and in the toy model of the last section, or are negligible, then the conclu- 
sion would be that the free-energy density of QED decreases significantly 
at super high temperatures, suggesting a phase transition |]IE||. This high- 



temperaure phase of QED might then be analogous to the low-temperature 
phase of QCD: one might have bound states of the electrons, positrons and 
photons. Or, the high-temperature phase might be due to the formation of 
other structures, such as magnetic strings ||18||. Of course at the moment 



all this is speculative as one knows neither the sign nor magnitude of the 
non-perturbative corrections. 



9 Conclusion 

A method has been developed which enables one to obtain bounds on the 
full pertubation expression, S(\), of a physical quantity eventhough the only 
information available is its partial perturbation series 

N 

S N (\) = Y,fn* n - (59) 

n=0 

The first requirement is that not all the f n be of the same sign. Then so- 
lutions exist to the extremum condition (|TT]). If the first nontrivial coefficient 
f n {n > 0) is negative, then the solutions to ( pTTj) will be global or local min- 
ima, depending on the sign of f^. The approximants Sn as defined through 
(|T3| ) then form lower bounds to S(X) if the inequality (|i~6l) is satisfied for 
all N (and if the plausible assumptions leading to (|15|) are admitted) . The 
inequality ([16| ) must first be explicitly tested for the available terms of the 
series, then the arguments given in the Appendix indicate that the trend will 
continue. Therefore even with partial information as in floTf), one can deduce 
lower bounds to the exact value S(X). Additional bounds may be obtained 
by using the auxiliary series method described in Sec. (5). 

Sometimes one finds minima solutions to (|Tl"D but for which the Sn obey 
an inequality opposite to (P^[). In such cases, even though the convergence is 



monotonic, it is not a priori obvious that the Sn are actually upper bounds to 
the exact value. For problems where the exact result is unknown, additional 
input from theory or physics is required before a definitive statement can 
be made. However in all the examples encountered of this type, the Sn did 
actually bound the exact result. A similar caveat concerns the approximants 
S defined near the end of Sec. (2). 
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The observed rapid convergence of the bounds has been explained in the 
Appendix. The obvious question is whether the bounds converge to the exact 
value itself. Empirically, it is found that the bounds formed by Sn actually 
converge to the exact value if S(X) has a slowly varying first derivative, 
dS(X)/d\. Otherwise the complementary bounds formed by 5(A) or the 
auxiliary series method give better approximations to the exact result. The 
reason for this phenomena has been given in the Appendix. 

Currently the most popular resummation methods for divergent series are 
the Pade' or Borel-Pade', see for example |I| [19| . Those methods usually do 
not converge monotonically, and so do not provide bounds, but nevertheless 
can sometimes be used to estimate the sum of a series. Those estimates can 
probably be constrained by using bounds obtained through the method of 
this paper. 

Although most of the observed features of the resummation method devel- 
oped here have been explained in the Appendix, at least in a semi-quantitative 
way, more patterns were detected than could be explained. In order to high- 
light some of these apparently universal trends, I summarize them as three 
questions: (i)Is it true that in all cases, as iV — > oo, c(N) = p(N+l)/p(N) — > 
1 (Similarly for the p{N)) ? (ii) Is it true that in all cases the approximants 
Sn and Sn form bounds to the full perturbative result ? (iii) Is it true that 
the bounds Sn always approximate the exact result S well if | /|| | is small 
throughout the range, < A < A , of interest, and otherwise the alternative 
bounds Sn or those from the auxiliary series give better approximations ? 

The technique itself can be improved in several ways. For example, if the 
large order behaviour of /„ is ~ (2n)!, instead of n\, then a generalised Borel 
transform can be used to take advantage of that fact. Alas, in order to keep 
this paper itself from growing out of bound, some of these refinements and 
further physical applications have to be discussed elsewhere [|T1J] . 
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Appendix 
Al 

The resummed perturbation series with a variable parameter p is given 

by 



n=N f f A\ n N ~ n (2n + k - IV r°° 

S »^^M Ziz&hwL ^ ( ,A)<->. (60) 



n=o nl \Pj t'o k\(2n-l)\ Jo 
Write the equation above in the compact form 



MM = £^, (6i) 

n=0 V 

where A n (N) is a p and A dependent constant that can be read off from flBOp . 
In particular note that the sign of A n (N) is the same as the sign of f n . The 
extremum condition flTT|) applied to fl6"T|) results in 

dS N » nA n (N) 

where use has been made of the form w(z\) in transforming the p derivative 
of A n (N) into a A derivative. At the solution p = p(N) of ([TlD , one has 

Since for a given N there is in general more than one solution to the 
extremum equation (|ll]) ? here p(N + 1) and p(N) will refer to solutions at 
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consequtive orders which are both positions of minima or both positions of 
maxima. Now define 

p(N + 1) = c(N) p(N) (64) 

where c(N) is some function of N. Then from equation ( ]6"2"| ) and the corre- 
sponding one at order N + 1, one easily deduces 

d(S N+1 -S N ) » _jx_ ( A n (N + l) _ \ (N + 1)A N+1 (N + 1) 

dX Uo ~[ p(N) n \ c{N) n n{ ') c{N) N + x p{N) N + 1 

(65) 

Now for large N, A n (N) ~ A n (N + 1), and write this simply as A n . Also 
define, ASn = SV+i — SV Assume now that p(N) is large and c(N) > 1. 
Also assume that all the coefficients A n are generically of the same order, or 
at least do not increase rapidly with n (the factorial growth of f n has already 
been taken care of by the Borel transform). Then keeping terms needed to 
solve for l/c(N) at leading order, ( |65|) simplifies to 



OAS* _A_ (J_ \ 2A^ 
d\ |Ao p(N) \c{N) J p(N) 2 ' { } 

Similarly, from eq.(|63|) and the corresponding one at order iV + 1 one deduces 
at large iV and to leading order in l/p(N), 

""-mi^)- 1 )- (67) 

where it is implicit in this equation that A = Ao- Now, if Ao is varied, then 
the leading change in fl6"7|) comes from A\. Thus using (0), Eq.(p6]) can be 
simplified to 

a ( 1 ,\ 2A 2 



"I (68) 



p(iv) V c ( Ar ) / p( n ) 

where a = (A\ — A^-)|a - Thus, 

_L ~ i + 2 ^ 2 (69) 
c(N) + ap(N) ' 1 J 

This solution will be self-consistent with the initial assumption c(N) > 1 
only if A 2 /a is negative, which then requires that that f 2 and f\ should be of 
opposite signs (since A\ is larger than j^)- Note that for N = 2 the solution 
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p{2) exists in the first place only if f 2 and fi are of opposite signs. Hence 
one concludes that c(2) > 1 is generically expected. Comparing ( |69|) with 
the corresponding equation at the next order iV + 1 one obtains at large N 
the recursive relation 

1 11 

+ , (70) 



c(N + l) c{N) c{N) 

which can also be written as 

1 1 



c(N+l) c(N) V °( N ), 



(71) 



showing that 

c(N + 1) < c(N) . (72) 
Indeed the large N solution of (|70| ) is 

' 1 + Tr~~~r? ■ (73) 



c(N) N + K 

with K a constant. Therefore c(N — > oo) — > 1 + , that is, though p(N) will 
increase with N, it will approach a constant as iV — > oo. Now, since p(N) 
increases with N, this means that the various approximations that led from 
to ffffltffflj will become increasingly accurate. What is remarkable is that 



in the examples studied with c(N) > 1, ( [72]) is already satisfied at low 
and furthermore the relation (|7TJ|) too is a reasonable approximation. 

The above relations (|70| - |72D were obtained under the assumption c(N) > 
1. If c(N) < 1 then there are apparently no simple relations. For the example 
in Sec. (3.5), one had c(N) < 1 for all N. For the auxiliary series of the beta 
function in Sec. (7), c(N) alternated between being slightly larger than one 
and much smaller than one. 

Consider now the Eq.(^) which is valid at large iV and large p(N), 

^-jfefe- 1 )" (74) 

Remarkably, this simple equation summarizes most of the observed trends. 
When c(N) > 1, (|74|) implies that for A% < 0, which corresponds to f\ < 
and minima solutions to ( PD , ASn > 0, which is indeed observed in the 
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examples studied and leads to lower bounds. Similarly if A\ > 0, which cor- 
responds to f\ > and maxima solutions, ASn < 0, implying upper bounds. 
The only exception observed so far is the Euler-Heisenberg series in Sec. (6) 
which had c(N) > 1, A\ < 0, and yet gave upper bounds. Presumably for 
that example the approximation (|74]) is not appropriate. Now suppose that 
c(N) < 1, as happens in the model of Sec. (3.5) and (roughly) for the auxil- 
iary series in Sec. (7). Then Eq.([74|) implies that minima, corresponding to 
fi < and hence Ai < 0, give AS < 0, which explains some of the oddities 
observed. 

The equation ( |74|) also explains the rapid convergence of the bounds. 
Indeed one sees that convergence can be achieved simply by one of two ways. 
Firstly, if c(N) — > 1 as N — > oo, that is p(N) — > po- This type of behaviour, 
is manifested, for example, by the toy-model in Sec. (3.5). The second way 
is for c(N) > 1 for all N, which leads to p(N) increasing with N. This is 
what has been generically observed. Actually, as shown above, the second 
behaviour also gives c(N) — - > 1 + , and hence the convergence is doubly fast. 



Explicitly, from Eqns.( |69| , [73| , |7^ ) one deduces for the c(N) > 1 case, 



AS ~ ^L(J--l) 2 (75) 
2A 2 \c(N) J 1 ; 

2A 2 (N + K) 2 ' 1 ' 

In summary, for practical problems, the strongest statements can be made 
for the case when c(N) > 1 is observed for the given terms of a partial 
perturbative series, that is when p(N) increases with N. Then the various 
approximations leading to the above equations become increasingly accurate 
at large N, allowing one to make assertions about all N. Firstly, one deduces 
c(N — > oo) — > 1 + . Secondly, for Sn which are global (and local) minima, 
if Sn < Stv+i for the given terms of the series, the trend will continue for 
larger N, the convergence of the Sn will be rapid, and they will form lower 
bounds to the exact result. However if Sn > SV+i is observed for the global 
(and local) minima, (as in Sec. (6)), then though the trend will continue and 
though the convergence of the Sn will be rapid, it is not a priori obvious 
that they will be upper bounds to the exact result, because then key pieces 
(|TJ-|15D of the argument are missing. (The same loophole occurs for the 
approximants Sn formed from the p(N)'s.) 

If it is observed that c(N) < 1 for the given terms of a series, then the 
various equations above are not necessarily accurate at higher N. In that 
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case one can only conjecture that the observed monotonic and rapid conver- 
gence will continue at higher orders. 



A2 

Consider now the slope of the approximants Sjv(A). From (|60|), as A — ► 0, 
Sat(A) -> /o + /iA, so that 

fW/, (TT, 

Thus all the bounds approach the origin with the same value (/o) and slope 
(/i) independent of p and iV.This fact can be seen in all the figures. Consider 
next Eq.QB2{) for large p(N), 

^U„J_J**_ (78) 
d\ 1X0 X oP (N) ' K } 

Since the sign of A 1 is the same as the sign of fi, this shows that the curves 
Sjv(A) are not expected to change direction as A varies. This is indeed ob- 
served, and explains why the bounds Sat (A) to S(X) are also good estimates 
of S(X) itself only when the latter has a slowly varying first derivative. 

As was observed in the main text, the approximants S(X) formed from 
the local extrema p(N) had a more varying slope. In order to understand 
this, set for simplicity fo = so that Sjv(O) = and let us demand that 

S N (X ) = 0, (79) 

so that Sat (A) curves back to its value at the origin, and so can better ap- 
proximate functions like sin(A7r). The condition ([79]) then leads from fl6l|) 
to 

N A 

= E^ (80) 
^1 P( N ) 

Note that since the p(N) also have to satisfy the extremum condition (|TTD, 
eqns. ( p^pO|) are actually two coupled equations for Ao and p(N) which we 
would like to analyse for consistency. Firstly, (|S0|) can be used to eliminate 
the leading term in the 'barred' version of (|62|), so that now 

9S„ _ " (n - l)A n 
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which at large p(N) gives the slope 

dS N , 1 A 2 



d\ lXo \ p(N)< 



(82) 



If /2 is opposite in sign to f% then by comparing (J82D with ( |77D one sees that 
indeed the approximants S can change direction Of course this just shows 
that the approximate large p(N) analysis above is self-consistent. Now, in 
Sec. (2), for fi < 0, f%> 0, the p(N) have been defined as positions of local 
maxima when the p(N) are positions of local minima. We can compare the 
relative magnitudes of the two values as follows. For large p(N), (|8*0|) has 
the approximate solution 

P(N) ~ ^ ■ (83) 

By contrast the p(N) are approximate solutions of (p2]) with the left-hand- 
side deleted, which is equivalent to ignoring the mild p dependence of the 

—9 A 

p(N) ~ . (84) 

Thus the values of p(N) are expected to be smaller than those of p(N), and 
the interested reader may verify from the given examples that this is indeed 
the case. Reversing the logic of the argument above, one concludes as follows. 
For that iV when p(N) is the position of a local minimum, one expects a local 
maximum at p(N). While the Sn curve is expected to be monotonic in A, 
that of 5(A) will not be monotonic if the value of p(N) is smaller than that 
of p(N). 

A similar discussion can be carried out for the curves formed from an 
auxiliary series S N . Again set for simplicity f = 0, and define 

& - ^ • (85) 

The extremization fllTl) in p is done with respect to the resummed auxiliary 
series S N , and one deduces for large p (N), from an equation analogous to 



(H), that 



Now demanding 



sign^A^] =sign(/ 2 ). (86) 



S N (X )=0, (87) 
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for the reconstructed resummed series Sn = XSj^ gives, 




( 



as: 



'N 



) 



(88) 



Ao - 



which when combined with ( |36"D shows that for fi and f2 of opposite signs, 
the slope of the reconstructed Sn at A = Ao is opposite in sign to its slope 
at the origin which is given by flTTp. 

Thus if f\ and / 2 are of opposite signs, then the auxiliary series may be 
expected to give a reconstructed Sn with a slope that has variation in sign, 
compared with the slope of the approximant Sn which is obtained by direct 
means, if the two conditions (|TTD and (|87D for S have a consistent solution 
p'(N) for some Aq. 
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Figure Captions 

Figure (la): Plots of the divergent series Sn(X) for the model in Sec. (3.1), 
together with the exact result S(X). Starting from the lowest curve and mov- 
ing upwards, one has S5, S3, S, S2, S4. 

Figure (lb): Plots of the resummed series Sjv( A) for the model in Sec. (3.1), 
together with the exact result S(X). Starting from the lowest curve and mov- 
ing upwards, one has S2, S3, S 4 , S. 

Figure (lc): Plots of the resummed series S^ N (X) for the model in 
Sec. (3.1), together with the exact result S(X). Starting from the lowest 
curve and moving upwards, one has £4, S, S 5 , S 3 . Compared to Fig. (lb), 
the convergence is faster but not monotonic. 

Figure (2a): Plots of the divergent series Sn(X) for the model in Sec. (3. 2), 
together with the exact result S(X). Starting from the lowest curve and mov- 
ing upwards, one has S 6 , S 2 , S, S4. 

Figure (2b): Plots of the resummed series Sn(X) for the model in Sec. (3. 2), 
together with the exact result S(X). Starting from the lowest curve and mov- 
ing upwards, one has S4, S 5 , S 6 , S 7 , S 8 , S. 

Figure (3a): Plots of the series Sjv(A) for the model in Sec. (3.3), together 
with the exact result S(X). Starting from the lowest curve and moving up- 
wards, one has S3, S, S4, S2. 

Figure (3b): Plots of the resummed series Sjv(A) for the model in Sec. (3. 3), 
together with the exact result S(A). Starting from the lowest curve and mov- 
ing upwards, one has S2, S3, S4, S. 

Figure (4a): Plots of the divergent series Sat (A) for the model in Sec. (3.4), 
together with the exact result S(A). Starting from the lowest curve and mov- 
ing upwards, one has S 3 , Si, S, S 4 . 

Figure (4b): Plots of the resummed series Sjv(A) for the model in Sec. (3.4), 
together with the exact result S(A). Starting from the lowest curve and mov- 
ing upwards, one has S 4 , S 5 , S 6 , S 7 , S. 
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Figure (5a): Plots of the divergent series Sjv(A) for the model in Sec. (3. 5), 
together with the exact result S(X). Starting from the lowest curve and mov- 
ing upwards, one has S3, S, S 2 , S4. 

Figure (5b): Plots of the resummed series Sat (A) for the model in Sec. (3.5), 
together with the exact result S(X). Starting from the lowest curve and mov- 
ing upwards, one has S, Sq, S4, S 2 . The approximants approach the exact 
result from above. (The curves for Sq and S are indistinguishable). 

Figure (6a): Plots of the series SV(A) for the model in Sec. (3.6). Starting 
from the lowest curve and moving upwards, one has Si, S 5 , Sg, S 7 , S3. 

Figure (6b): Plots of the resummed series Sjv(A) for the model in Sec. (3. 6), 
together with the exact result S(X). Starting from the lowest curve and mov- 
ing upwards, one has S3, S4, S5, Sq, Sj, S. 

Figure (6c): Plots of the series Sjv(A) for the model in Sec. (3.6). Starting 
from the lowest curve and moving upwards, one has S, Sq, S5. The approx- 
imants approach the exact result from above. 

Figure (7a): Plots of the divergent perturbation series S N (X) for the 
model in Sec. (4). Starting from the lowest curve and moving upwards, one 
has S 5 , S 3 , S 2 , S 4 . 

Figure (7b): Plots of the resummed perturbation series Sat (A) for the 
model in Sec. (4), together with the full perturbative result S pert (X). Start- 
ing from the lowest curve and moving upwards, one has S 2 , S 3 , S 4 , S 5 , S pert . 

Figure (7c): Plots of S exact (X), Ss(A) and Ss(A) — S np (A) as defined in 
Sec. (4). At A = 10, the lowest curve is of S5 — S np , and the highest one is S5. 

Figure (7d): Plots of the resummed perturbative result S 2 (A) and the 
exact (sum of perturbative and non-perturbative) value S exact (X) as defined 
in Sec. (4). At A = 2, the lower curve is of S 2 . 

Figure (8): Plots of the resummed series Sat(A) for the model in Sec. (5.1), 
together with the exact result S(A). Starting from the lowest curve and mov- 
ing upwards, one has S, S 7 , Sq, S 5 . The approximants approach the exact 
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result from above. 

Figure (9): Plots of the resummed series <Sjv(A) corresponding to the aux- 
iliary series of (pE8| ) in Sec. (5.2), together with the exact result S(X). Starting 
from the lowest curve and moving upwards, one has S, Sq, S5, S4, 5*3, S2. 
The approximants approach the exact result from above. 

Figure (10): Plots of the resummed series Stv(A) corresponding to the 
auxiliary series of (|50|) in Sec. (5.2), together with the exact result S(X). 
Starting from the lowest curve and moving upwards, one has S, Sq, S$, S4. 
The approximants approach the exact result from above. 

Figure (11a): Plots of the divergent Euler-Heisenberg series Sjv(A) given 
in Sec. (6). Starting from the lowest curve and moving upwards, one has 

1S3, £>2> S4. 

Figure (lib): Plots of the resummed Euler-Heisenberg series Sjy(A), to- 
gether with Schwinger's exact result S(X). Starting from the lowest curve 
and moving upwards, one has S, Sq, S4, S2. The approximants approach 
the exact result from above. 

Figure (12a): Plots of the resummed beta function iSjv(A) of Sec. (7). 
Starting from the lowest curve and moving upwards, one has Sy, Sg, 5*5, 04, S3, 5*2. 

Figure (12b): Plots of the resummed series Sjv(A) for beta function 
in Sec. (7). Starting from the lowest curve and moving upwards, one has 
S3, S 5 , S 7 . The approximants appear to form upper bounds. 

Figure (12c): Plots of the resummed beta function <Sjv(A) °f Sec. (7), 
obtained through the auxiliary series. Starting from the lowest curve and 
moving upwards, one has S2, S3, S4, S5, Sq. The approximants appear to 
form upper bounds. The curves for N = 2 and N = 3 are indistinguishable, 
and similarly, those for N = 4 and N = 5 are very close. 

Figure (12d): Magnification of the N = 6 curve of Fig. (12c) near its non- 
trivial zero. 

Figure (13a): Plots of the divergent perturbative free-energy density of 
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QED, Sjv(A) given in Sec. (8). Starting from the lowest curve and moving 
upwards, one has S 2 , S 5 , S 3 , §4. 

Figure (13b): Plots of the resummed perturbative free-energy density of 
QED, Sat (A). Starting from the lowest curve and moving upwards, one has 
5*3) S4, S5. 
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Fig. 7b 
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